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Abstract 

We compute the gluino lifetime and branching ratios in Split Supersymmetry. Using 
an effective-theory approach, we resum the large logarithmic corrections controlled by 
the strong gauge coupling and the top Yukawa coupling. We find that the resummation 
of the radiative corrections has a sizeable numerical impact on the gluino decay width 
and branching ratios. Finally, we discuss the gluino decays into gravitino, relevant in 
models with direct mediation of supersymmetry breaking. 



1 Introduction 



The long gluino lifetime is a trademark of Split Supersymmetry PJ El IS] • The experimental 
discovery of a slowly-decaying gluino [4J would not only be a strong indication for Split 
Supersymmetry, but it would also allow for a measurement of the effective supersymmetry- 
breaking scale m, which cannot be directly extracted from particle dynamics at the LHC. 
Moreover, the gluino lifetime is a crucial parameter to determine the cosmological constraints 
on the theory |TJ|3]. Therefore, for both experimental and theoretical considerations, it is 
very important to have a precise prediction of the gluino lifetime and branching ratios. 

For what concerns the gluino decay processes in the MSSM, tree-level results for the 
decays into chargino or neutralino and two quarks and one-loop results for the radiative 
decay into neutralino and gluon can be found in the literature jHj. In Split Supersymmetry, 
however, the quantum corrections to the gluino decay processes can be very significant, 
because they are enhanced by the potentially large logarithm of the ratio between the gluino 
mass rrig and the scale m at which the interactions responsible for gluino decay are mediated. 
A fixed-order calculation of these processes in Split Supersymmetry would miss terms that 
are enhanced by higher powers of the large logarithm. In order to get a reliable prediction 
for the gluino decay width, the large logarithmic corrections have to be resummed by means 
of standard renormalization group techniques. 

Recently, a calculation of the gluino decay widths in Split Supersymmetry was presented 
in ref. working at tree level for 3-body decays and in (not resummed) one-loop approx- 
imation for 2-body decays. In this paper we will present a calculation of the gluino decay 
processes that includes all leading corrections in a s and a t , the strong and top- Yukawa cou- 
pling constants. As we will show, the inclusion and resummation of leading-order corrections 
give sizeable modifications of the gluino branching ratios, even for moderate values of m. 

The structure of the paper is as follows: in sect. El we list the operators in the effective 
Lagrangian of Split Supersymmetry that are responsible for the decays of the gluino, and 
the high-energy boundary conditions on the corresponding Wilson coefficients; in sect. OJwe 
determine the renormalization group evolution of the Wilson coefficients, and we express 
the operators in the low-energy effective Lagrangian in terms of mass eigenstates; in sect.EI 
we discuss our numerical results for the branching ratios and total width of the gluino 
decays in Split Supersymmetry; in sect. El we consider the possibility of gluinos decaying 
into gravitino; in sect.|H]we present our conclusions. Finally, in the appendix we provide the 
analytical formulae for the gluino decay widths. 
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2 The Effective Lagrangian 



Below the squark and slepton mass scale to, the effective Lagrangian of Split Supersymme- 
try describes the dynamics of Standard Model (SM) particles together with higgsinos and 
gauginos. At the level of renormalizable interactions, there is a conserved G-parity (under 
which only the gluino is odd) preventing gluino decay. However, integrating the squarks out 
of the underlying supersymmetric theory induces non-renormalizable interactions that vio- 
late the G-parity. Restricting our analysis up to dimension-6 operators, the G-odd effective 
Lagrangian at the matching scale to is given by 

c = ^j:c°Q* + ^- 2 jzc?Qt + 4 (j2 c ?Q? + h - c ) • (i) 

to 2 ^ to ~^ m 2 J 

We are working in the basis of interaction eigenstates for gauginos and higgsinos, ne- 
glecting the effect of electroweak symmetry breaking, since to ^> M z . The G-odd operators 
involving the B-ino (B) are 



Qf = BYl 5 ~9 a ® E q2'h l *T a qW (2) 

k=i 

Ql = B^ l5 ~g a ® E 7(i r«W (3) 

fc=i 

Qf = BYl 5 9 a ® E dM^VdW (4) 

fc=l 

Q? = ^7^75^ ® ?i 3) 7^ a ?i 3) (5) 

Q 5 5 = l7"75^° ® t Rl »T a t R (6) 

Q 6 5 = l7"7 5 ^ ® Kl^T a b R (7) 

Q* = I<7^ 75 5°^, (8) 

where k is a generation index, T a are the SU(3) generators and G^ v is the gluon field strength. 
Assuming that the squark mass matrices are flavour-diagonal, the Wilson coefficients of the 
operators Qf at the matching scale to are 

Cf(m) = Cf(m) = - 9 -^f r u , Cf{m) = C}{m) = ^ r fijJ , (9) 
C 3 5 (m) = CH^) = r dfl , G 7 5 (to) = ^ (to, - to~) E (r& - r fa ) Q s , (10) 

where = to 2 /to~. Note that G 7 B vanishes for mass-degenerate squarks. 
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The G-odd operators involving the WMno (W) are 
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QY = W A Yl5 9 a ® E q ( L k) I 1 * r A T a (11) 



fc=i 



= W*-f lb g a ® # 7 M r A T a gf, (12) 
where r A are the Pauli matrices. The matching conditions for the Wilson coefficients are 

C?(m)=C?(m) = - 9 -fr, L . (13) 

For the higgsinos, we use a compact notation in which the two Weyl states H u and Hd 
are combined in a single Dirac fermion H = H u + e H d , where e is the antisymmetric matrix 
(with £12 = 1) acting on the SU(2) indices. The states H u and H d can be recovered by chiral 
decomposition, H u = H L and H d = — e {H C ) L . Keeping only the third-generation Yukawa 
couplings, the G-odd operators involving higgsinos are 

Qf = H L ? R ® eq^T a t R (14) 

Qf = H L a^~g a R ® eqfa^T a t R (15) 

Ql = Th~g a L ® b R T a q^ (16) 

Q} = J^ R a^~g a L ® b R a^T a q^ (17) 

Ql = H L a^g a R hG% (18) 

where h is the Higgs doublet. The Wilson coefficients at the matching scale m are 

C?{m) = (r & - r, R ) , G/(m) = |^ (r & + r fi J , (19) 

V2 sin /? 4 V 2 sin p 



C"{m) = -J^- (r SL - r dR ) , G/(m) = - % h ( r& + rjj , (20) 
V 2 cos p 4v2cosp 



Qf(™) = o 9 d h l ■ J r*. + (2i) 

32 v 2 7r sm p 

Here foj and /if, are the top and bottom Yukawa couplings, and tan (3 is a free parameter of 
Split Supersymmetry. 

Before proceeding to the operator renormalization, we want to make some remarks. 

(i) We recall that all coupling constants appearing in the expressions of the Wilson 
coefficients given above have to be computed at the scale m. 

(ii) Note that we have given the Wilson coefficients of the 4-fermion operators at the 
leading perturbative order, while the coefficients of the operators Qf and Q" are given 
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at the next order (one-loop approximation). The operator anomalous dimensions will be 
computed in sect. [3] at the leading order in the strong and top- Yukawa couplings a s = 
glKA-n) and a t = ft|/(47r). Therefore, the gluino 3-body decays, mediated only by 4-fermion 
operators, will be computed by resumming all a Sj jln(m/mg) corrections, but neglecting 
terms 0{a^ 1 ln n {m/rrig)] with n > 0. For the radiative 2-body gluino decay into a gluon 
and a neutralino, a greater accuracy is more appropriate. The expressions of Cf and C" 
given in eqs. (|1U|) and ([21]). together with leading-order anomalous dimensions and one-loop 
matrix elements [see eq. (|62j) below], allow us to determine the 2-body decay amplitude 
neglecting terms 0[a™^ 1 ln n (m/m§)] with n > 1. This means that we have resummed all 
large logarithms at the leading order in all cases, but our formulae for 2-body gluino decays 
contain also the complete 0(a s j) terms, relevant when the logarithm is not large. 

(Hi) If m is close to the GUT scale, in presence of gauge-coupling unification there is 
no solid justification for the approximation of computing a s contributions to the anomalous 
dimensions, neglecting electroweak corrections. However, because of the large SU(3) coeffi- 
cients, we consider our approximation to be fairly adequate, even for m as large as 10 13 GeV, 
which is the maximum value of m consistent with the negative searches for anomalous heavy 
isotopes. 

(iv) In eq. fl2())) we have included the contribution from the bottom Yukawa coupling hb, 
since these coefficients are enhanced when tan {3 is large. There are no tan f3 enhancements 
in the evolution below m, and therefore our results are reliable for any value of tan/3. 

(v) Split Supersymmetry is free from flavour problems, therefore our assumption that 
squark mass matrices are diagonal is unnecessary. On the other hand, a certain degree of 
mass degeneracy among squarks is required by gauge-coupling unification. In the results 
presented in sect. EJwe take for simplicity all squark masses to be equal. 

3 Operator Renormalization 

The renormalization-group flow for the Wilson coefficients is determined by the equations 




2tt ' 



(22) 



(23) 



(24) 
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where /i is the renormalization scale and, in Split Supersymmetry, we have (3 S = 5, (3 t = —9/2 
and j3 st = 8. The anomalous-dimension matrix 7 can be expressed as 



nrij = -2bi j -S ij ^2af, (25) 
/ 

where bij are extracted from the poles of the one-loop renormalization of the operators Qi 
(Qi — > Qj/e + • • •). In eq. (|23|) the sum is over all fields entering the operator Q iy and the 
field anomalous dimensions a/ are given by 

a q k = -j— (oi s C F + 7^3) , a u k = (a s C F + a t 5 k3 ) , a dR = - ^7— - , (26) 

L 4 7T V Z / fl 4 7T 47T 

a ^ = -^T' ah = -^T> a * = ^{ Nc -3 N 0- (27) 

Here /c is a generation index, Cp = (iV? — 1)/(2N C ), N c = 3, Nf = 6. Note that the gluon 
anomalous dimension a g (given here in the Feynman gauge) is different from the SM value 
because it includes the gluino contribution. 

We find that the anomalous-dimension matrices of the 5-ino operators in eqs. (|2*|)— (|H)l . 
of the W-ino operators in eqs. (fiT|) - (fr2*|) . and of the higgsino operators in eqs. (|T4^ - (|THj) are 

respectively 

= 4vr s 4vr * 4?r 7^ a = b,w,h (28) 
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(32) 



(33) 



For coefficients with only multiplicative renormalization (which is the case for Cf , G\ 
C^ 3 4 ), eq. (|2*2^) can be easily integrated, with the result 



w 
2- 



Ci(ii) = Ci(m) r]. 
We have defined 



/ 7s Pstlt ] 
\2f3 s 2/3 a ft J 



Vt 



It 



for d = C 7 B , C{ 



2- 



/n H 
° 1,3,4- 



Vt 



a s {m) 

Oitim) 
a t (/i) 



a s {m) fi 

1 + — ^ Ps In — , 

2tc m 



a t {m)(3 t 



r]s s + 



a s (m) (fist - (3 S ) 



(34) 

(35) 
(36) 



The evolution of the Wilson coefficients for the other 5-ino operators involves operator 
mixing and the solution of eq. (f2*2*|) is given by 



cm 

cf(») 



rfs 



_9_ r 

10 



Cfim) + yC(m)\ i = 1,2,3,6 
\l + z) Cf(m) - z Of (m) + y C{m) 
[1 + 2z) Cfim) - 2z Cfim) + yC{m 



(37) 
(38) 
(39) 

where C = C?/3 + (Of + C* + C 4 i )/6 + (C£ + C 6 5 )/12, y = r/ s 4 / 5 - 1, and z = {v! /15 Vt 1/3 ~ 
l)/3. Because of the non-vanishing contribution from r j" tl the equations for C" and C" 
cannot be solved analytically The numerical results for the renormalization coefficients A^-, 
defined by 

(C*Q*)\ (A 22 A 25 \(C«(m)\ 

\C*M) U 52 A 5 J{ C »(m)) ' 1 Uj 

are shown in fig. for a representative choice of a s (m) and at(m). Despite the fact that the 
high-energy boundary condition on C", eq. (jUJ, is suppressed by a loop factor, a sizeable 
value of C 5 H can be generated through the mixing with C 2 . 
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Figure 1: Renormalization group flow of C" and C 5 H , expressed in terms of the coefficients Ay 
of eq. ()40|) . for a s (m) = 0.05 and a t (m) = 0.03. The solid, dashed, dotted, and dot-dashed 
lines correspond to A 2 2, A 2 5, A 52 and A 55 , respectively. 

A computation of the 0(a s ) part of the anomalous dimensions, restricted to the four- 
fermion operators, has been given in the appendix of ref. [Sj. From the comparison with 
eq. (f2HJ) it appears that the authors of ref. [Sj have omitted the mixing among the 5-ino 
operators induced by the penguin diagrams. Also, we disagree with ref. jS] on the anomalous 
dimensions of the higgsino operators. 

Once we have evolved the Wilson coefficients down to the renormalization scale at which 
we compute the gluino decay width, it is convenient to express the operators in terms of 
chargino and neutralino mass eigenstates. With the usual definitions for the chargino and 
neutralino mixing matrices U, V and N, which we assume to be real, the .B-ino, jy-ino and 
higgsino spinors can be expressed as 

W+ = xt {Un Pl + V iX P R ), H^ = xf (U i2 P L + V l2 P R ) , (41) 

B = ^N iU W*=^N l2 , T°=xl(NuPR-N i3 P L ), (42) 

where Pl and Pr are the chiral projectors. In the basis of mass eigenstates, the effective 
Lagrangian becomes 

£ = =iE p/ ? Q? + =2 °f of + h - c -) • ( 43 ' 
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The operators involving neutralinos and quarks and their corresponding Wilson coeffi- 
cients are 

2 



k=l 



QsLvn = ^h,l91 r ® q R , L T a q L , R 
QaLo r = ^ R ,^ V l^l R ®q R ^^T a q L 



(q 


= u,d) , 


(44) 


(q 


= t, 6) , 


(45) 


(q 


= t, 6) , 


(46) 


(q 


= *, 6) , 


(47) 





= CfN a 


+ C*N a , 




--C 2 B N a , 




v° 
^ld L 


= CfN a 







--C°N a , 




v° 


= C*N tl 


+ C?N l2 , 




/^fXj _ 
°3t L - 




X° H 

C 4 t L = C 2 Ni4 , 


U 2t R 


= C*N tl 






r x» _ 

^4t R - 


—C 2 N i4 , 


v° 


= C»N iX 






C«N i3 , 


C 4 b L = -C 4 iVj 3 


U 2b R 


= C e fiV ll 




_ 


-C~ H N lZ 


v° 

C<M — 
°46 fl - 





(48) 
(49) 
(50) 
(51) 
(52) 
(53) 



The operators involving charginos and quarks and their Wilson coefficients are 

QiL = tf., R l*r L , R ®i:d? lfl T a uP (54) 

k=l 

Q2L = xt L , R r~g a L , R ® ^7j a tL (55) 

Q 3 X L = xf^fiE,* ® b R , L T a t L , R (56) 

<?4t* = ld R ^ u ~9tn®K L ^T a t LiR (57) 



Cft = -V^C^i, = V2C*U iU (58) 

C 2 X 2 = — y/2 C 2 Va , C3I =C"Ui 2l C*l =C 4 H Ui 2 , (59) 
C 2 x l = V2 Cf U a , C*l = Cf V i2 , Cfl = Of V l2 . (60) 

All Wilson coefficients in eqs. (|4*Hj) - (|5^j) and (|S8*|) - (|5n|l are evaluated at the scale // at which 
the gluino decay width is computed (we take \i = nig in our numerical analysis). 

The magnetic operator involving a neutralino and a gluon is 

=^^1^0^. (61) 

In order to reach the desired accuracy in the g — > gx° process, we need to include the matrix 
element contribution coming from the diagram in which the two top quarks in the operator 
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Q 2 H close in a loop emitting a gluon. This results into an "effective" Wilson coefficient 

° 9 X L (A*) = C?(aO N a + Of M ^ + f^T <3?(a0 ^4 v In ^ , (62) 

where t> is the Higgs vacuum expectation value and we take \i = m g . 

From the effective Lagrangian of eq. (j4*Hj) we can compute the gluino decay widths and 
complete expressions can be found in the appendix. The same effective Lagrangian correctly 
describes also the interactions that lead to the decays g —>■ ggx° and g — > gh°x°. However, 
since these processes are subleading, we will not explicitly calculate their decay widths. 



4 Results 



We are now ready to discuss the results of our computation of the decay width and branching 
ratios of the gluino in Split Supersymmetry. The input parameters relevant to our analysis 
are the sfermion mass scale m, the physical gluino mass rrig and tan/3, which in Split Su- 
persymmetry is interpreted as the tangent of the angle that rotates the finely tuned Higgs 
doublets. To simplify the analysis we assume that the squark masses are degenerate, i.e. we 
set rq L = Tu R = = 1 in the matching conditions of the Wilson coefficients. The gluino 
mass parameter in the Lagrangian, M 3 , is extracted from m g including radiative corrections, 
and the other gaugino masses Mi and M2 are computed from M3 assuming unification at 
the GUT scale. The higgsino mass parameter \x is determined as a function of M2 by re- 
quiring that the relic abundance of neutralinos is equal to the dark-matter density preferred 
by WMAP data [H] (see fig. 11 of ref. [2] ). The sign of fi remains a free parameter, but 
since it does not affect our results for the gluino decays in a significant way we will assume 
11 > throughout our analysis. The effective couplings of gauginos and higgsinos at the weak 
scale, needed to compute the chargino and neutralino mass matrices, are determined from 
their high-energy (supersymmetric) boundary values by means of the renormalization-group 
equations of Split Supersymmetry, given in ref. j2j. Finally, the SM input parameters rele- 
vant to our analysis are: the physical masses for the top quark and gauge bosons, m t = 178 
GeV, Mz = 91.187 GeV and Mw = 80.41 GeV; the running bottom mass computed at the 
scale of the top mass, mb{m t ) = 2.75 GeV; the Fermi constant, Gp = 1.166 x 1CT 5 GeV -2 ; 
the running strong coupling computed at the scale of the top mass, a s (m t ) = 0.106. 

To start our discussion, we show in fig. Elthe gluino lifetime Tg (in seconds) as a function 
of the sfermion mass scale m, for tan (3 = 2 and four different values of the physical gluino 
mass (rrig = 0.5, 1, 2 and 5 TeV, respectively). It can be seen that Tg is about 4 seconds 
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m (GeV) 



Figure 2: Gluino lifetime r g as a function of the sfermion mass scale m, for different values of 
the physical gluino mass m g . The other free parameters are chosen as tan/3 = 2 and \i > 0. 
The dashed horizontal line corresponds to the age of the universe, t v = 14 Gyr. 

for nig = 1 TeV and m = 10 9 GeV. A value of r g equal to the age of the universe (14 Gyr) 
corresponds to m = (1.1, 2.1, 4.5, 13) x 10 13 GeV for m g = 0.5, 1, 2 and 5 TeV, respectively. 

In fig. El we show the branching ratios for the three decay processes g — > g — > x QQ 
and g — > x ±( 1Q' (summed over all neutralino or chargino states) as a function of m, for 
tan/3 = 20 and three different values of m g : 500 GeV (upper plots), 1 TeV (middle plots) 
and 2 TeV (lower plots). The value of tan/? has little impact on these results. The plots on 
the left of fig. El represent the results of our full calculation, including the resummation of 
the leading logarithmic corrections controlled by a s and at- The plots on the right represent 
instead the lowest-order results that do not include the resummation. We obtain the latter 
results by replacing the Wilson coefficients of the four-fermion operators in the low-energy 
effective Lagrangian with their tree-level expressions in terms of gauge and Yukawa couplings 
[eqs. (fn])-([Tf)j). (fT3^) and (fTTj|) - (f2T)|) ]. and the Wilson coefficient of the magnetic operator with 
its one-loop expression. The plots in fig. El show that the branching ratio of the radiative 
decay g — > x°9 decreases for increasing m g and increases for increasing m. In fact, as stressed 
in ref. the ratio between the two-body and three-body decay rates computed at lowest 
order scales like mf/m"- [1 — \n(rn 2 /mf)] 2 , where the logarithmic term comes from the top- 
stop loop that generates the magnetic gluino-gluon-higgsino interaction. For large values of 
m, the resummation of the logarithms becomes necessary. Comparing the plots on the left 
and right sides of fig. El we see that the resummation of the leading logarithmic corrections 
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Figure 3: Branching ratios for the gluino decay channels x°9 (dashed lines), x°QQ (dotted) and 
X ±c L<i l (dot-dashed), summed over all possible neutralino or chargino states, as a function of m, 
for three values of mf 500 GeV (upper plots), 1 TeV (middle) and 2 TeV (lower). The curves in 
the left (right) plots do (do not) include the resummation of the leading logarithmic corrections. 
Other parameters are tan/? = 20 and \i > 0. 
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Figure 4: Effect of the radiative corrections on the partial widths for the decays g — > x°9 (dashed 
lines), g — * (dotted) and g — > x ±( 1Q' (dot-dashed) as a function of m. The parameters are 
chosen as m„ = 1 TeV, tan/3 = 20 and [i > 0. 



tends to enhance the three-body decays and suppress the radiative decay. The effect of 
the corrections on the branching ratios is particularly visible when, like in the middle and 
lower plots, neither the two-body nor the three-body channels are obviously dominant in 
the range 10 s GeV <m< 10 13 GeV, relevant to Split Supersymmetry. 

To further illustrate the effect of the resummation of the leading logarithmic corrections, 
we plot in fig. 0] the ratio T/Tq of the partial decay widths with and without resummation, 
for the processes g — ► x°9> 9 ~ > X°QQ an d 9 — > X^ QQ '■ We fix rrig = 1 TeV, tan (3 = 20 and 
/x > 0, but we have checked that the qualitative behaviour of the corrections is independent 
of the precise choice of the parameters. It can be seen from fig. HJthat for large enough values 
of m the radiative corrections can be of the order of 50-100%, and that they enhance the 
widths for the three-body decays and suppress the width for the radiative decay. 

To conclude this section, we discuss the scaling behaviour of the gluino lifetime and total 
decay width. The lifetime Tg = h/T to t can be written as 



Tn 



4 sec 



m 



ITeV 



x 



10 9 GeV 



(63) 



where the normalization iV is of order unity and depends on m and m g (and only very mildly 
on tan (3). In fig. El we show N as a function of m for tan f3 = 20 and three different values of 
the physical gluino mass {m g = 0.5, 1 and 2 TeV, respectively). The non-vanishing slope of 
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Figure 5: The normalization N of eq. ()63|) as a function of the sfermion mass scale m, with 
(solid lines) and without (dashed lines) resummation of the leading logarithmic corrections. The 
upper, middle and lower sets of curves correspond to m g = 0.5, 1 and 2 TeV, respectively. The 
other free parameters are chosen as tan/3 = 20 and fi > 0. 

N represents the deviation of the total gluino decay width from the naive scaling behaviour 
r tot oc m| / m 4 . The solid lines in the plot represent the results of our full calculation, whereas 
the dashed lines represent the lowest-order results that do not include the resummation. 
For low values of m g the contribution of the radiative decay dominates (see fig. EI); thus the 
total decay width departs visibly from the naive scaling and is significantly suppressed by 
the resummation of the radiative corrections. On the other hand, for large values of m g the 
three-body decays dominate, and the effect of the resummation is to enhance the total decay 
width. Finally, for the intermediate value m g = 1 TeV there is a compensation between the 
corrections to the radiative decay width and those to the three-body decay widths, and the 
net effect on the total decay width of the resummation of the leading logarithmic corrections 
is rather small. 



5 Gluino Decay into Gravitinos 

Split Supersymmetry opens up the possibility of direct tree-level mediation of the original 
supersymmetry breaking to the SM superfields, without the need of a hidden sector j^j. 
In usual low-energy supersymmetry, this possibility is impracticable: for F-term breaking 
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some scalars remain lighter than the SM matter fermions, and for D-term breaking gaugino 
masses cannot be generated at the same order of scalar masses. In Split Supersymmetry a 
large hierarchy between scalar and gaugino masses is acceptable, and indeed models have 
been proposed IB] with direct mediation of D-term supersymmetry breaking. 

Therefore, in Split Supersymmetry the original scale of supersymmetry breaking y/F, 
which is related to the gravitino mass by 

ft* F 

m3 / 2 = VT^V (64) 

could be as low as the squark mass scale m. This means that the interactions between the 
gluino and (the spin-1/2 component of) the gravitino, which are suppressed by 1/F, could 
be as strong as those considered in the previous sections, which are suppressed by 1/m 2 . 

For 777.3/2 <C mg, the gravitino interactions can be obtained, through the supersymmetric 
analogue of the equivalence theorem ^H] , from the goldstino derivative coupling to the super- 
current. This approximation is valid as long as \f~F~ -C 6 x [m g /ITeV) 1 / 2 x 10 10 GeV. Using 
the equations of motion, we can write the effective goldstino (G) interactions for on-shell 
particles as 

C = - \-m\ L q L q L - m\ R q R q R + g* a^ l5 G a \ G + h.c. (65) 

Below m, the effective Lagrangian describing the interactions between the gluino and the 
goldstino becomes 

^=4i;c?g?, (66) 
b i=i 

Qi = GYl59 a ® E q {k) ^T a qW (67) 

k=l,2 
q—u,d 

Ql = G^7 5 9 a ®qi 3) ^T a qW (68) 

Qi = GYi 5 g a ® t Rl ,T a t R (69) 

Q? = G^ l5 g a ® b Rl ,T a b R (70) 

Qi = Ga^ l5 g a G%. (71) 

The Wilson coefficients at the matching scale m are 

°i — °2 — °3 — °4 — — ' 5 ~~ ~~ 2 \[2 



Note that the coefficients of the interactions in eq. (jbbj) have no dependence on m, because the 
squark mass square in the propagators of the particles we integrate out is exactly cancelled 
by the squark mass square in the goldstino coupling in eq. (fBHjl. 



14 



The operator renormalization proceeds analogously to the discussion in sect. El The 



anomalous dimension matrix of the operators in eqs. (J57j) - (J7T|) is given by eq 



with 



Is 



/ IK _ q/v 


16 




16 


16 


4 


4 — 9iV 




4 


4 


2 


2 


2 


-9N C 


2 


2 


2 




2 


2 - 9iV c 


V o 















/0 





0^ 






1 - 


-2 







if = 


-1 


2 





, tI = 




Vo o 





o) 












(73) 



2iV7 - 18NJ 



0. 



(74) 



The evolution of the Wilson coefficients for the goldstino operators has the simple analytic 
form 



T] s 



_9_ 

10 



C % G {m)+yC G {m) 



1,4 



(i + z)c 2 G (fi)- 2 c 3 G H + i/C G (fi) 



r]7 z C 5 G (m) 



2^)C 3 G (m) 



2 2 C 2 G (m) +yC (m) 



8/15-1/3 



where C = 2^/3 + C 2 G /6 + (C 3 ° + C 4 G )/12, y = ^ - 1, and z = (r/f/ 15 r/ t " 
The quantities 7/ s and 774 have been defined in eqs. (133|) and (jHEJ), respectively. 



(75) 
(76) 
(77) 
(78) 
l)/3. 



The formulae for the gluino decay widths into goldstino and quarks and into goldstino 
and gluon can be found in the appendix. In fig.|Ulwe show the branching ratio for the process 
g — > G g as a function of the ratio y/F /fn, for m = 10 9 GeV, tan/3 = 2 , /i > and different 
values of the gluino mass. The branching ratio for the decay into goldstino and quarks, 
suppressed by phase space, is always at or below the 1% level. It can be seen from fig. |B1 
that the gluino decay into goldstino and gluon is largely dominant when \/~F is as low as 
m. In fact, the decays into charginos or neutralinos and quarks (relevant for large values 
of rrig) are suppressed by phase space, while the radiative decay into gluon and neutralinos 
(relevant for smaller values of rrig) is suppressed by mf/m'~ and a loop factor. With respect 
to the scaling behaviour outlined in eq. (|6*3*|) . the additional contribution to the total gluino 
decay width coming from the decay into goldstino and gluon can significantly suppress the 
gluino lifetime. In fact, for \fF = m the normalization iV in eq. (|6"3*|) takes on values of order 
40-50 for m > 10 8 GeV. 

On the other hand, the widths for the gluino decays into goldstino are suppressed by a 
factor m A /F 2 with respect to those for decays into charginos or neutralinos. Fig. |H1 shows 



15 






= 0.5 TeV 





= 1 TeV 


. _ m 

g 


= 2 TeV 


m 


= 5 TeV 


g 






Figure 6: Branching ratio for the decay g — > G g as a function of yF /m, for different values of 
the physical gluino mass m~ g . The other free parameters are chosen as m = 10 9 GeV, tan/3 = 2 
and \i > 0. 



that as soon as we depart from the condition yF = m the branching ratio for g — > G g falls 
off very quickly, and already for y~F jm as large as 10 the gluino decays into goldstino are 
below the 1% level. 



6 Conclusions 

If Split Supersymmetry is the correct theory to describe physics beyond the Standard Model, 
one of its most spectacular manifestations might be the discovery of a very long-lived gluino 
at the LHC In this paper we provided a precise determination of the gluino lifetime and 
branching ratios. Applying to Split Supersymmetry the effective Lagrangian and renormal- 
ization group techniques, we discussed the proper treatment of the radiative corrections that 
are enhanced by the large logarithm of the ratio between the sfermion mass scale and the 
gluino mass. We computed the anomalous dimensions of the operators relevant to the gluino 
decay, that allow us to resum to all orders in the perturbative expansion the leading loga- 
rithmic corrections controlled by a s and at- We also provided explicit analytical formulae 
for the gluino decay widths in terms of the Wilson coefficients of the effective Lagrangian of 
Split Supersymmetry. For representative values of the input parameters, we discussed the 
numerical impact of the radiative corrections and found that they can modify substantially 
the gluino decay width and branching ratios. Finally, we considered models with direct 
mediation of supersymmetry breaking, and we found that the gluino decays into gravitinos 
might dominate over the other decay modes. 
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Appendix 



We present in this appendix the explicit formulae for the leading three-body and two-body 
gluino decay widths. All the results are expressed in terms of the Wilson coefficients of the 
effective Lagrangian of Split Supersymmetry, discussed in sects. |2J El and El 



Three— body decays into quarks and chargino or neutralino: denoting the momenta 
of the decay products as {pi,P2,P3) = {Pq I -,Pq J ,Px)i an d Sij = (Pi+Pj) 2 , the three-body decay 
amplitude is given by 

1 



XQlQj 



256 7r 3 m~ m 4 



\M\ ds 13 ds 2 3 



(79) 



The bar over \A4\ denotes the average over colour and spin of the gluino and the sum over 
colour and spin of the final state particles (the dependence on m has been factored out). 
The limits of the integration in the (§13, S23) plane are 

1 



'13 



'13 



,max 
'23 

„min 



2 2"^ 

m n + m % H 

" * 2s 23 



4 



S23) («23 - "t|j + m 2 x ) 



(trig — m 

\ 1/2/ 2 2 \ \ 1/2/ 2 

+ A 1 (s 23 ,m~ g ,m qi )\ 1 (s 23 ,m^,m 

1 2 2 \ / 2 , 2 \ 

(m- - m qi - s 23 ) (s 23 - m lj + m x ) 

\ 1/2/ 2 2 \ \ 1/2/ 2 

-X' (s 2 3, rrig, m qi ) A (s 2 3,«V r , m 



2^ 

<?j' ""X^ 



qiJ j 



where A(x, ?/, z) 



(K 
(|m x | +m 5j ) 2 , 

x 2 + y 2 + z 2 — 2 [xy + xz + yz). 



(80) 

(81) 
(82) 
(83) 



In the computation of the decays involving quarks of the first and second generation we 
can neglect the quark masses and we find 



y y _ 

\fdu~ x~ud 



y _ 

X?<39 



1536 7r 3 m 4 



C& + C& \g(x i )-2CtiCftf(x i ) 



»4) 



v° 2 v° 2 

°l<?r "I" W 



- 3i? 



5f(xi) + f(xi) 



(q = u,d), (85) 



7687r 3 m 4 

where Xi = m Xi /m g , and we have included an overall factor 2 to take into account the sum 
over the two generations of light quarks. The functions / and g are defined as: 



g(x) = 1 ~8x 2 + 8x 6 -x 8 - 12 x 4 lnx 2 , 

f{x) = 2x + 18x 3 -18x 5 -2x 7 + 12a; 3 (l + a; 2 ) lnx 2 



(86) 
(87) 
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For generic quark masses the integration of the squared amplitude \M\ 2 on the (s± 3 , s 23 ) 
plane cannot be performed analytically, and in order to compute the total decay width we 
must resort to a numerical integration. 

The squared amplitude for the processes g — > xt bt and g — > x7 ^ i s given by 



|A4| 2 = (m? + m\ - s 13 ) (si 3 - m 2 + - m 2 ) 



+ C^h ( m s + m b - s 2s) («23 - m 2 + - m 2 ) 



+ 2 +2 

+ 7 ( C 3 X 2 + C 3 X A ) (™!+ + m 5 " S 13 " S 2 3 ) (Si 3 + S23 " ™t ~ 



+ 4 C 4 X < + C 4 X A 



(m 2 + + m~ - si 3 - s 23 ) (si3 + s 2 3 -m 2 - ml - 4m 2 + ) 



+4 (si 3 - m 2 + ) (S23 - ™ 2 + ) - 4 m 2 m 2 

A-i A. j 



+ 2 C 2 X 2 C2 X A m § m x + ( s i3 + S23 - m 2 y + - m~) 



+ ( C 2 X i C 3 X i + 12 C 2 x i Cll ) m + m, (s 23 - m\ - m 2 ) 



+ [Ct Cll C* m § m b (s 23 -m 2 -m 



- Ct C& - 12 Ct Ct R m~ g m t (s 13 - m 2 b - m\ 



- (C 2 X / C 3 X / - 12 C x / C x /) m x + m b (s 13 - m 2 - m 2 ) 



+ 2 c& c# +c 3 ^a 



IR 



(m~ + m 2 + - si 3 - s 23 ) (s 23 - S13 + m 2 - m 2 ) 



1 o 2 / 2 2 \ <-> 2 / 2 2 \ 

+2 m 6 (s 23 - m t - m +) -2m t (s 13 - m b - m +) 



2 ( C& + 48 C 4 X 2 C X A ) m 5 m t m b . 



(88) 
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The squared amplitude for the processes g — > "X^tt and g — > x° & & is given by 



v° 2 v° 2 



'2<ta 



(m? + m 2 - s 13 ) (s 13 - ml - m 2 ) 



+ (m? + m 2 - s 23 ) (s 23 -ml- m 2 ) + 2m g mo (m~ + m 2 - s 13 - s 23 ) 



1 /^ x o 2 2 



(m Y o + m~ - s 13 - s 23 )(si 3 + s 23 - 2 m 2 , 



, Y° 2 V° 2 



m x° + m l + 4m 9 _ s i3 - «23)(si3 + s 23 - 2 m 2 , - 4m 2 o' 



+4 (si 3 - m 2 , - m 2 o) (s 23 - m 2 - m 2 o) + 8 m 2 m 2 o) 



"9 ""Xi 



+ 4 (7 2 Xl L (7 2 x ^ m 2 (s 13 + s 23 + 4 m 5 m x o - 2 m 2 ) 



+ 



v o v o v o v o 



m x o (m? + m 2 - s 13 ) + m~ g (m 2 o + m 2 - s 23 ) 



+ 



L '2g J? ^3^ _ °2g L ^3<7 fl j m q 



m x o (ml + m 2 - s 23 ) + m g (m 2 + m 2 - s 13 ) 



m„ 



m- g (m 2 o + m 2 q - s 23 ) - m x o (m| + m\ - s 13 ) 



x° x° x° x° 
12 ( C 2 g R C 4 g L - C 2 ^ C^ R ) m q 



m x o (ml + m 2 - s 13 ) - m g (m 2 + m 2 - s 23 ) 



2 f C 3 x l Cl\ L + C$ B Cl\ R ) (m? + m 2 + 2 m 2 . - s 13 - s 23 ) (s 23 - s 13 ) 



2 (^ L c , 3 x ( ; fl + 48^ r a Xi 



(q = t,b). 



19) 



We have checked that inserting in eqs. (|HHj) - (jH^j) the high-energy (i.e. non resummed) 
expressions for the Wilson coefficients given in sects. |2] and El we reproduce the tree-level 
results of ref. HI. 
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Two— body decays into neutralino and gluon: the width for the radiative decay of the 
gluino, g -> g x - , is 

(m? -m 2 ) 3 2 

r n = g * f cl*? ) . (90) 

The use of the effective coefficient C s g ff defined in eq. (jo^jl allows us to reproduce the 
complete one-loop result when the resummation is switched off. 



Decays into goldstino: the decay width into goldstino and quarks of the first and second 
generation is: 

77? — 2 

r~ 9 C? , (91) 



G io 192 vr 3 F 2 
where we have summed over all four light quark flavours. 

The gluino decay width into goldstino and third-generation quarks is as in eq. ()79j) . with 
m replaced by F 2 . The squared decay amplitude, which has to be integrated numerically 
on the (S13, s 23 ) plane, is given by 



\M\ 2 = [C°-+C°-] (m 2 - + m 2 -srs)(sn-m 2 )+ (m 2 + m 2 -s 23 )(s 23 -m 2 q 



2 ~ 2 



+ 4C q G L C q G R m 2 q (s 13 + s 23 -2m 2 q ) (q = t,b), (92) 

where 

Cf h =cl=ch C? R = C°, Cl = cf. (93) 



Finally, the gluino decay width into gluon and goldstino is: 



Gg 2 7tF 2 



^ Cf . (94) 
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